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We develop a theory of Coulomb oscillations in superconducting devices in the limit of small 
charging energy Ec <C A. We consider a small superconducting grain of finite capacity connected 
to two superconducting leads by nearly ballistic single-channel quantum point contacts. The tem- 
perature is supposed to be very low, so there are no single-particle excitations on the grain. Then 
the behavior of the system may be described as quantum mechanics of the superconducting phase on 
, the island. The Josephson energy as a function of this phase has two minima which become degen- 

^\ • erate at the phase difference on the leads equal to tt, the tunneling amplitude between them being 

0^ I controlled by the gate voltage at the grain. We find the Josephson current and its low-frequency 

■ fluctuations and predict their periodic dependence on the induced charge Qx = CVg with period 2e. 
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Q ■ I. INTRODUCTION. 

in ■ 

Coulomb effects in several different types of tfiree-terminal devices consisting of an island connected to external 
^-H . leads by two weak-link contacts, and capacitatively coupled to an additional gate potential, have been extensively 
^ studied during last years. The systems with a normal- metal island and leads were studied theoretically both in the 
tunnel-junction limit ||l| and in the case of a quantum point contact with almost perfect transmission The theory 
' of charge-parity effects and Coulomb modulation of the Josephson current was investigated in details in ^ . All the 
, above-mentioned systems at present are realized experimentally. Recently it was shown to be possible to produce 
T-H quantum point contact between two superconductors via a normal-conductive region made of two-dimensional electron 
. gas (2DEG) g; smeared step -wise behaviour of the critical current was observed, in qualitative agreement with 

■ predictions Q for the superconductive quantum contact with a few conduction channels of high transmittivity. An 
I observation of a non-sinusoidal current-phase relation in superconducting mechanically-controllable break junctions 

■ has been reported in |^ , again in agreement with . Another interesting experimental achivement was reported in , 
^ , where S-N-S contact with a size comparable to the de Broghle wavelength in the N region made of BiPb was realized 

'■^ ' and nonmonotonic behaviour of the critical current with the thickness of normal region was found. This remarkable 
Ch , development of technology points to the principal possibility to make a system of a small superconductive (SC) 
island connected to the superconductive leads by two quantum point contacts (QPC). In such a system macroscopic 
quantum effects due to competition between Josephson coupling energy and Coulomb (charging) energy could be 
realized together with quantization (due to small number of conductive channels) of the Josephson critical current. 
In the present paper we develop a theory for an extreme case of such a system, namely, for the case of two almost 
^ ballistic one-channel QPCs connecting a small SC island with two SC leads. We consider the limit of the characteristic 
charging energy much smaller than the superconducting gap, Eq <SC A, and, therefore, the Coulomb effects are small. 
We derive the dependences of the average Josephson current across the sytem, and its fluctuations (noise power) as 
functions of the SC phase difference between the leads a, and of the electric gate potential Vg. The Coulomb effects 
reveal themselves at phase differences a close to tt, when the two lowest states are almost degenerate. We show that 
such a system realizes a tunable quantum two-level system (pseudo-spin 1/2) which may be useful for the realization 
of quantum computers (see e.g. ). 

The paper is organized as follows. We start with considering a single QPC connecting a superconducting island 
to a single lead (Section II). We find the oscillations of the effective capacitance on the island as a function of the 
gate potential (in some analogy with Matveev's results |^ for a normal QPC). Depending on the backscattering 
probability in the contact, it may be described either in adiabatic or in diabatic approximation. We find the condition 
for the diabatic-adiabatic crossover. Then in Section III we formulate a simple model for the double-contact system 
in the adiabatic approximation. We replace the full many-body problem by a quantum-mechanical problem for the 
dynamics of the SC phase on the middle island. In Sec. IV we calculate average Josephson current through the system 
as a function of a and Vg, with a particular emphasis on the case of the phase difference a close to tt (when our 
effective two- level system is almost degenerate). Sec.V is devoted to the analysis of the Josephson current noise; we 
calculate integrated intensity So of the "zero" -frequency noise (an analogue of the noise calculated in [p2|-p^ for a 
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single superconduvtive QPC) as well as finite- frequency noise S^^ due to transitions between the two almost-degenerate 
levels. Finally, we present our conclusions in Sec. VI. 

II. ADIABATIC-DIABATIC CROSSOVER IN A SINGLE SUPERCONDUCTING QUANTUM POINT 

CONTACT. 

Consider a small superconducting island connected to an external superconducting lead by an one-channel nearly 
ballistic quantum point contact PJT^]. The electric potential of the grain may be adjusted via a gate terminal (fig. la). 
Following 1^ we assume that the contact is much wider than the Fermi wavelength (so that the transport through 
the constriction may be treated adiabatically) , but much smaller than the coherence length = Jivf/ttA (where vp 
is the Fermi velocity, A is the superconducting gap). 




FIG. 1. (a) Single QPC. The system consists of a SC grain connected to a SC lead via a QPC. A gate terminal is used to 
control the electric potential of the grain, (b) Double-contact S-S-S system. The second terminal is added to the single-QPC 
setup. 

Our assumption of low temperature is that the average number of one-electron excitations on the island is much 
less than one. Then they cannot contribute to the total charge of the grain and we may restrict our Coulomb 
blockade problem to the evolution of the superconducting phase only. The condition of low temperature is then 
T < A/ log(yz^(0)A), where V is the volume of the grain, j/(0) is the density of electron states at the Fermi level. 

We neglect phase fluctuations in the bulk of the island and describe the whole island by a single superconducting 
phase X- At a fixed value of the phase on the island, the spectrum of the jimction consists of the two Andreev states 
localized on the junction and the continuum spectrum above the gap A jl^] (fig. 2). The energies of the Andreev 
states lie below the gap: 

E{x) = ±A^l-tsin\xm, (1) 
where x is the phase difference at the contact, t is the transmission coefficient. 



E 




FIG. 2. Single-contact energy spectrum. The spectrum consists of the continuum of delocalized states and the two Andreev 
(subgap) states. Dashed lines denote Andreev states in the absence of backscattering (diabatic terms). Solid lines are the states 
split by backscattering (adiabatic terms). 

At t = 1, the spectrum of Andreev states (|l|) has a level crossing point at x = t^. At this point, the left and right 
Andreev states have equal energies, but in the absence of backscattering {t ~ \) the transitions between them are 
impossible. Therefore, we expect that an ideal ballistic contact cannot adiabatically follow the ground state as the 
phase X changes, but remains on the same left or right Andreev state as it passes the level-crossing point x = tt. We 
borrow the terminology from the theory of atomic collisions and call the (crossing) Andreev levels at t = 1 diabatic 
terms (dashed lines in fig. 2), and the split levels — adiabatic terms (solid lines in fig. 2). Instead of transmission 
coefficient i, it will be more convenient to speak of the reflection coefficient r = 1 — t. At r = 0, the contact is 
described by diabatic terms. As r increases, the transitions occur between the terms, and at sufficiently large r the 
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system will mostly adiabatically follow the split Andreev levels. In this section we study adiabatic-diabatic crossover 
and find the crossover scale for the reflection coefficient r. 

We assume that the reflection probability r <C 1 (almost unity transmission) and that the charging energy Ec <C A 
(the charging energy is defined by Ec = (2e)^/C). The latter assumption appears natural, because, like in the tunnel 
junctions | p7[ | we expect that the capacitance C of the grain has an additional contribution from the capacitance of 
the point contact. This capacitance is of order A/e^. A more detailed discussion of this phenomenon will be given 
elsewhere. At the moment we just mention that this contribution to the capacitance leads to the inequality Eq < A. 

To probe the degree of adiabaticity, we study the periodic dependence of the ground state energy Eq on the gate 
voltage. Because of the weakness of charging effects, this dependence will be sinusoidal: 

EoiVg) = ecos{2TrN) (2) 

(where N — VgC/2e is the dimensionless voltage), and we are interested in the amplitude e of these oscillations. The 
physical meaning of this periodicity is the oscillations of the induced charge on the grain — it follows immediately 
from the relation 

There is a simple physical explanation of the sinusoidal dependence (||). The ground-state energy modulation is 
determined by phase-slip processes in the contact. Such processes are phase tunneling events with phase changing 
by ±27r. While the magnitudes of the clockwise and counter-clockwise tunneling amplitudes are the same, their 
phases are ±27rA^. This results in the expression (^. Higher-order tunneling processes would give rise to higher-order 
harmonics in the periodic iV-dependence. This argument shows that the amplitude of oscillations e coincides with 
the phase-tunneling amplitude and, therefore, provides a good measure of adiabaticity in the phase dynamics. 

Under assumption Eq <C A, we may describe the contact by the dynamics of the phase on the grain and thus 
reduce the problem to a single-particle quantum mechanics. Since we restrict our attention to low lying excitations, 
it is only necessary to include the two Andreev levels on the junction. The potential term is the Josephson energy of 
the Andreev levels, the kinetic term is the charging energy. After a simple computation of the backscattering matrix 
elements (the off-diagonal entries in the potential term) , we arrive to the following Hamiltonian: 



H = H{x) + -EciTT^ - iV)2 (4) 



where 



Here x is the phase difference across the contact, r is the reflection coefficient. Obviously, the eigenvalues of H{x) 
reproduce the result (^. The number of Cooper pairs at the grain tt^ is the momentum conjugate to Xi [Xi'^'xl ~ *• 
Notice that x takes values on the circle x = X + 27r, and, accordingly, tt^ is quantized to take integer values. We may 
also write tt^ = —id/dx- 

This Hamiltonian loses its validity at the top of the upper band at x = 27rn, where the upper Andreev state mixes 
with the continuous spectrum (fig. 2). Howerver, the probability of the phase x to reach the top of the upper band 
of H{x) is exponentially small at Ec ^ A (smaller than the tunneling probability). The adiabatic-diabatic crossover 
is determined by the properties of the system near the minimal-gap point x — Therefore, we may neglect the 
transitions to continuous spectrum at x = 27rn. At the same time, we must disregard tunneling porcesses via the top 
of the upper Andrees band (ncxt-ncarcst-neighbor tunneling) which is present in the Hamiltonian (^-(^), but not 
in the original system. The nearest-neighbor tunneling is a feature of our model and is beyond the precision of our 
approximation. 

There are two opposite limits of the problem: small and "large" reflection. 

At zero reflection, the Hamiltonian splits into lower and upper components. Within each component the potential 
is periodic with the period An. As explained above, we must neglect the next-nearest-neighbor tunneling via the top 
of the bands. Therefore, the potential minima of H{x) are disconnected and cannot tunnel to each other, e = 0. 

The opposite limit is the case of "large" reflection (the precise meaning of " large reflection" consistent with r <C 1 
will be clarified below). In this limit, the gap opens in the spectrum of Andreev states, and the system adiabatically 
follows the lower state. We can replace the two-level Hamiltonian H{x) by its lowest eigenvalue and arrive to the 
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quantum-mechanical problem of a particle in a periodic potential. The quasiclassical limit of this problem is solved 
in the textbook [[l8|. In our notation the answer reads as follows: 



Ead = const V-Ec A exp ( - S'c/ ) , 



where 



Sc 



Bi 




— log 

4. ^Ec 



0(1) 



(6) 



(7) 



is the classical action connecting two nearest minima (or more precisely the two return points). The numerical 
constant Bi is of order one (at r 0, Bi — 4.69 + 1.41rlogr + ...). 

To study how the adiabaticity is destroyed it is useful to introduce the dimensionless "coherence factor" /(r) defined 

by 



fir)£ad, 



(8) 



where Sad is the amplitude of oscillations of the ground-state energy derived in the adiabatic approximation (with 
only the lowest Andreev state included). We see that /(O) — 0, f{r ^ Vad) = 1. The crossover scale rad can be 
derived by computing the corrections to /(r) in these two limits. 

First consider the limit of weak backscattering (r <^ rad)- In this limit we take the wavefunction to be the ground 
state of the Hamiltonian with zero r (at a given wavevector N), and then compute the first-order correction in r^^^ 
to the energy. The wavefunction is of "tight-binding" type and is generated by the "ground-state" wavefunctions ^'i 
localized in the potential minima (diabatic terms). The components of the two-dimensional vectors alternate: 



*.(x) 





i+1 







Then we find 



2(*,|^fi2(x)l*.+i) -2ri/2A I dx**(x)*.+i(x)sin 



(9) 



(10) 



(We assume the wavefunctions to be normalized). It is important that ^'^ and are wavefunctions for different 
potentials (— Aq cos(x/2) and Aq cos(x/2)) and the overlap integral (10) has a saddle point at the minimal-gap point 
X = TT, and it reduces the effective region of integration to Ix"'"'! ^ {Ec /^Y^^ ■ The normalization of the quasiclassical 
tail of the wavefunctions ^'i(x) yields 



*(x = 7r) =exp(-S',;(x = 7r)) 
(up to a numerical factor independent of Ec/A). Thus we obtain 



i.e., in terms of the "coherence factor" /(r), 



fir) 



r.l/2 



A_ 

E^ 



1/4 



(11) 



(12) 



(13) 



The physical meaning of the integral ([l0| ) is the summation over all paths shown in fig. 3a. 





FIG. 3. Tunneling paths in the diabatic (a) and adiabatic (b) limits. These diagrams represent the lowest-order corrections 
to the phase tunnehng amplitudes in the diabatic and adiabatic Umits respectively. 
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The above calculation shows that the crossover scale to adiabatic behavior is 



p \ 1/2 



rad-\^-f) . (14) 

In fact, we neglected the effect of change in the classical action Sd due to opening a gap; this effect is estimated to 
be of order 




(15) 

i. e. it is a higher-order effect than the change in /(r) proportional to r^/^. Notice that the characteristic scale of 
this change in the classical action is again Tad ^ yEc/^ (corresponding to 5Sci ^ 1). 

We may alternatively find the crossover scale rad by computing the lowest order correction to the "coherence factor" 
/(r) in the adiabatic limit. In this limit the Hamiltonian (^,^ may be rewritten in adiabatic terms (the voltage N is 
for simplicity moved to the boundary condition '^{x + Stt) = e^*'^^5'(x) by a gauge transformation) as 

H - -f (|)^ + Dix) _ f [G(x)| + |g(x)] - f G^(x), (16) 



where 



D{x)-['''^^^ A.) (17) 



E^ix) 



is the diagonalized form of the matrix (|^), and 



G(x)=(_;(^) 'o'^O' 

and |0) and |1) are the eigenvectors of the matrix (^. The last term in the Hamiltonian ( p^ ) can be shown to give 
smaller corrections than the term of the first order in G(x)- A careful perturbation theory in g{x) gives in second 
order 

l-f{r)^ [ e^ite,>^^)-^^te,'^^)5(Xi)5(X2)dXidX2, (19) 

where S'i,2(xijX2) ^^re the classical actions along the lower and the upper adiabatic branches between the points xi 
and X2- This integral corresponds to summation over all tunneling paths shown in fig. 3b. The function g{x) for the 
given matrix H[x) is a lorentzian peak at x = tt of height r~i/^ and width r^/^. Putting everything together, the 
integral ( p^ ) is calculated to be 

This asymptotics agrees with the found previously crossover scale (^) . 

To summarize the results of this section, the characteristic scale for adiabatic-diabatic crossover in a nearly-ballistic 
single contact is found to be rad ~ yEcl^- The phase tunneling amplitude is proportional to the gate-voltage 
modulation of the effective capacitance of the island, and thus can be directly measured. At low reflection coefficients, 
these oscillations are proportional to -/F, like in the normal 1-channel QPC B. 



III. ADIABATIC APPROXIMATION OF A DOUBLE-JUNCTION SYSTEM. 

Now turn to the case of a double-junction system (fig. lb). As before, we assume that the reflection probabilities 
in both contacts are small, ^ 1, that the charging energy Ec <C A and that the temperature is sufficiently low 
to prohibit single-electron excitations on the grain. To adjust electrostatic potential of the grain we again use a gate 
terminal, N = VgC/2e denotes the dimensionless gate voltage, as before. 

For the moment, to simplify the discussion we assume that the reflection coefficients in the contacts are greater 
than the crossover scale rad found in the previous section and, therefore, we may consider only the lower adiabatic 
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branch of the Andreev states. In fact, the results may be extended further to the case < rad by using appropriate 
"coherence factors" /(r), similar to those in the previous section. 

We set the superconducting phase on one of the leads to be zero; the phase on the other lead a is assumed to be 
fixed externally. Then the total Josephson energy of the two contacts is (fig. 4): 

U{x)^Ui{x) + U2[a-x), (21) 

where 

U^{5<t)) = -Ay^l-tj sin2((5</>/2) (22) 
are the lower adiabatic Andreev terms in the two junctions. 



ViX) 








FIG. 4. Potential U{x)- At a 7^ it has two minima 
potential, but leaves the bottom of the wells unchanged. 



a/2 K+al2 1 

Finite backscattering in the contacts smoothes the summits of the 



At ti = ^2 = 1, the potential U{x) obviously has two minima — at % = q;/2 and at x = + tt — and sharp 
summits at x = cind % = tt + a (fig. 4). At small nonzero r,, gaps open at the crossing points of Andreev levels, 
which smoothes the summits of U{x)- Still, the bottom of the potential remains practically unchanged. 

The adiabatic Hamiltonian for the double junction looks like follows: 

H{a, N) - U{x) + U{a - x) + \ec{-i^ - Nf. (23) 

2 dx 

The potential term of the Hamiltonian is the sum of Josephson energies of the contacts, the kinetic term is the 
Coulomb energy of the charge at the grain. 



IV. JOSEPHSON CURRENT. 



The condition Ec ^ A allows us to treat the Coulomb term in the Hamiltonian perturbatively. First, neglecting 
the Coulomb term, we obtain a classical system on the circle in the potential ( |2T| ) with two minima. The energies of 
the minima are Vi{a) = — 2A| cos(a/4)| and V2{a) = — 2A| sin(a/4)| (see fig. 4). To a very good precision, we may 
neglect backscattering in determining the minima — except near the point a = 0. Since all the Coulomb effects occur 
near the resonance point a = tt, this approximation is justified. At zero temperature, our classical system prefers the 
lowest of the minima. Thus the energy of the S-S-S system in the absence of the Coulomb term is given by 

£;(«) = -2Acos(a/4) for - vr < a < tt (24) 

(see fig. 5). Differentiating this energy with respect to the phase a gives the Josephson current 

I (a) = 2e ^^}"^ = A sin - for - tt < a < tt (25) 
da 4 

(fig. 6). Notice that the current has large jumps at the points of level crossing a = tt + 27rn. Qualitatively this picture 
is very similar to the case of a single S-S ballistic junction, but the shape of the current-phase dependence I{a) is 
different. 
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FIG. 5. Classical minimum of the potential U{x) as a function of the external phase difference a. Dotted line shows the 
quantum gap opened by the Coulomb term. 
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FIG. 6. Josephson current as a function of the external phase difference a. Dotted line shows smearing of the singularity 
due to the Coulomb term. 



If we assume a non-zero temperature T <C A, the occupation of the upper minimum is exponentially small except 
in the vicinity of the level-crossing point ja — 7r| ~ T/A. Thus, the effect of the temperature results in the smearing 
of the singularity in I{a) ai a — tt. 

Another source of level mixing near the singular point a = tt is quantum fluctuations, i.e. the fluctuations arising 
from the kinetic term in the Hamiltonian (^3|). They result in nonzero amplitudes of tunneling through the two 
potential barriers between the potential minima. Due to the shift in the "angular momentum" by N, the wave 
functions in the two potential wells aquire an additional factor exp{iNx)- This results in the relative phase of the 
two tunneling amplitudes by 2t:N. The net tunneling amplitude (defining the level splitting) may be written as 

HMN) = A7(iV) = A(7ie"^ + j^'^''). (26) 

where 71 and 72 are the two amplitudes of phase tunneling in the two different directions (i.e. of phase slip processes 
in the two different contacts). Below we assume that these amplitudes are computed at the level-crossing point a = tt, 
where they are responsible for level splitting. 

The amplitudes 71 and 72 obey all asymptotics derived in the previous section (except for numerical factors) . When 
the backscattering in the contacts r ^ rad, they may be found in the quasiclassical approximation: 




71,2 ~ -T eM-B2\ -E^) « 1, (27) 



where i?2 ~ 1 is determined by the classical action connecting the two potential minima (at r ^ 1, B2 = 1.45 -I- 
2.20rlogr + ...). At r <C Vad, the tunneling amplitudes are 

7i,2'-ri/2exp(-B2\/^) (28) 

For the best observation of Coulomb oscillations, 71 and 72 must be of the same order, but not very small. In the 
ideal case 71 = 72 = 7 the total amplitude 

7(iV) = 27COs(7riV) (29) 
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Although the periodic dependence ( p9| ) has 4e period as function of the "external charge" Qx = CVg = 2eN ^ the 
Josephson current and its fluctuations depend on |7(A^)p only (cf. eqs.( ^p4| ) below), and their period is 2e as 
expected ||]. 

The characteristic scale for the r-dependence of B2 is 5r ^ ^jEcjK^ therefore for 71 and 72 to be of the same 
order, the transparencies of the two contacts must differ by no more than \r\ — r2\< \J Ec / 

Here we should comment on the difference of our result (|26|)-(|2^) from the normal two-channel system discussed 
in H]. In the normal system the two tunneling amplitudes multiply, and the net ground-state energy oscillations are 
proportional to rlnr at small r. In the superconducting system, the external leads have different superconducting 
phases, and the tunneling in the two contacts occurs at different values of the phase on the grain. Therefore, the 
tunneling amplitudes add with some phase factors and give the asymptotic of ^Jr at r — > 0. In fact, the oscillations 
in the superconducting system will be proportional to r (similarly to the normal system ^) in a different limit — at 
the phase difference a — Q, when the potential U{x) has a single minimum and a single barrier. 

The hybridized energy levels in the vicinity of a = tt are given by the eigenvalues of the 2x2 Hamiltonian 



H{a,N) 

Diagonalization gives the two energy levels: 



Ei^2ia,N) = -A 



sm — N- cos — ± 
4 4 



. a 
|cos- 



72 (TV) 



(30) 



(31) 



the off-diagonal matrix elements of the Hamiltonian open a gap at the level-crossing point a — tt (fig. 5). This gap 
periodically depends on the gate voltage Vg, and these oscillations comprise the Coulomb effects in the S-S-S junction. 

We can obtain the Josephson current by differentiating the energy levels with respect to the phase a. The gap 
results in smearing the singularity in I (a) even at zero temperature (fig. 6): 



/(a) = ^sm(^^ 



V^sin2(^) + i72(iV) 



for a 



(32) 



The width of the crossover at a = tt depends periodically on ¥„: \a ■ 



In the above discussion we neglected the excited oscillator states. The interlevel spacing for the excitations in the 
potential wells is of order ^/ AEc S> A7. Therefore the Coulomb effects have a much smaller energy scale and the 
excited states do not participate in mixing the ground states of the two potential wells. 

At a nonzero temperature these Coulomb effects will compete with the smearing by temperature so that the width 
of the singularity at q = tt is given at nonzero temperature T <C A by |a — 7r| ~ max(7(A^), T/ A). Therefore, in order 
for Coulomb effects to dominate the thermal fluctuations, we must have T < 7A. 

It is instructive to compare this picture with the case of multi-channel tunnel S-S-S junction (to distinguish from 
the results of we should remark that we consider the opposite to their assumption A < Ec limit). If we develop a 
similar theory for tunnel Josephson junctions, we find that the potentials (^l|), ( p^ are both sinusoidal, and, therefore, 
the total potential ( |2l| ) has only one minimum (versus two in the nearly ballistic system). In the tunnel S-S-S system 
the current-phase relation I{a) has a smearing at a — n due to the difference between the critical currents of the 
two Josephson contacts. The Coulomb effects compete with this smearing and in order to win, the charging energy 
Ec must be greater than the difference of the critical currents. In the tunnel system the corresponding splitting 7 
is linear in Ec while in the nearly ballistic system it is exponentially small. Otherwise, Coulomb oscillations in I{a) 
will appear similar in these two cases. 

To summarize this section, we observed that the Coulomb effects in the one-channel S-S-S junction smears the 
singularity in the Josephson current I{a) at the critical value a = (2n + l)7r. This smearing depends periodically on 
the potential of the grain with the period 2e/C and is exponentially small in the adiabatic parameter Eel A ^ 1. 
The smearing is the result of mixing the two states in the potential minima of the Josephson energy. 



V. FLUCTUATIONS OF THE JOSEPHSON CURRENT. 



In this section we compute the low-frequency spectrum of the fiuctuations of the Josephson current in our model. We 
shall be interested in frequencies much less than the oscillator energy scale \fKEc^ thus we consider only transitions 
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between the eigenstates of the reduced ground-state Hamiltonian (|30|). We also assume that the temperature is lower 
than ^/AEc, then we may disregard the excited oscillator states and the internal noise in the contacts (discussed in 
|p^- p^ , p^ ) . Obviously, under these assumptions we can observe current fluctuations only in the close vicinity of the 
resonance point a = ±7r, where the energies (^l]) of the two low-lying states are close to each other. 

We expect to observe two peaks in the noise spectrum — one at zero frequency (due to the thermal excitations 
above the ground state), and the other at the transition frequency \Ei — -E2I (from the off-diagonal matrix elements 
of the current operator). In this section we compute the intergal weights of these peaks and postpone the discussion 
of their width (determined by dissipative processes) until elsewhere. 

Discuss first the zero-frequency peak. In our approximation it is just the thermal noise of a two-level system. In 
the vicinity of the resonance point a = n we can linearize the spectrum 14,2(0^) and make an approximation that one 
of the two states carries the current I{a, N), and the other —/(a, N). The spectral weight of the noise is then given 
by a simple formula: 



So{a,N, T) 



{If 



cosh^ 



(33) 



Substituting /(a, N) and i?i_2(Q!, N) from the previous section, we obtain the noise intensity near the resonance: 

2 

-2 



5'o(a,iV,r) = 



A2 



a — TT 

2V2 



cosh 



A // a-Tr V 
TV UV2 j 



72 (TV) 



(34) 



For the effect of the Coulomb interaction to be observable, the temperature must be smaller than the Coulomb gap: 
T < 7A. At constant T and iV, the noise decreases exponentially as a goes away from its critical value a = tt, and 
at a = TT the noise is suppressed in the interval |a — 7r| < 7(A^) (fig. 7). Interplay between these two factors results 
in the strong dependence of the peak value of the noise on the potential of the grain. The peak value of the noise 
maXtt 5(q!, A^, T) is plotted against N in fig. 8. Most favorable is the case of identical contacts, when 71 = 72 = 7 
and, therefore, 7(A^) = 27Cos(7rA^). In this case, when cos(7rA^) ^ Tf-fA (small gap limit) the noise takes its 
maximal value S « In the opposite limit of large gap (cos(7riV) 3> T/7A) the noise decreases exponentially: 



S 



■ exp(— 4 



A7I cos ttA^I > 



LA7I cos TT ATI ^^^y T 

levels on the grain with different electron numbers have equal energies. 



The noise has a sharp peak at the resonance point cosTriV = 0, where two 




ycos ttN 

FIG. 7. Zero frequency noise as a function of the phase a. It decays exponentially for a far from the resonance point a - 
At the very resonance point, the noise is suppressed, because both of the two states carry nearly zero Josephson current. 




T/yA 



FIG. 8. 
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Maximal value of the noise versus the potential of the grain. The period of the peaks corresponds to the period 2e 



of the induced charge Q = CVg. The width of the peaks depends on the capacity of the grain. 
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Now turn to the noise peak at the interlevel frequency ui = \Ei — E2\- Since now cu can be large compared to T, 
one needs to discern between different kinds of frequency-dependent correlation functions, which can be measured 
as a noise intensity in different experimental situations here we mean by noise the Fourrier spectrum of the 
time-symmetric current-current correlation function. In our approximation of a two-level system such a noise is 
temperature independent, and its weight is determined purely by the off-diagonal matrix element: 



S -i 



;i|/|2) (35) 



A straightforward computation for the Hamiltonian (^0|) and / = 2e{dH/da) gives (in the vicinity of a = tt): 

(1|/|2) = ^!lM(cos^ + sin^) (36) 



/A7(7V)\ 



2 



a — TT 



SUa,N) = A'^^^j cos'^. (37) 



and 



This result contrasts the corresponding noise intensity in the single quantum point contact (found in |12| , |l9| , [l3[ ) . In 
the single quantum point contact the correponding noise intensity S^j is temperature-dependent, because that system 
has four possible states (or, alternatively, two fermion levels). In the case of the double junction the system has only 
two states differing by the phase on the grain, and the quantum fluctuations 5*1^ become temperature-independent. 



VI. CONCLUSIONS 



We have developed a theory of Coulomb oscillations of the Josephson current and its noise power via the S-S-S 
system with nearly ballistic quantum point contacts. The period of Coulomb oscillations as function of the gate 
potential is Vg — 2e/C. These oscillations arise from the quasiclassical tunneling of the superconducting phase on 
the grain and are, therefore, exponentially small in Ec / A at Ec ^ A. In addition, we predict a crossover from 
adiabatic to diabatic tunneling at the backscattering probability Tad ~ \/ Ec j A. At backscattering below fad, the 
amplitude e of the Coulomb oscillations is proportional to the square root of the smallest (of the two contacts) 
reflection probability yj Train ■ This constrasts the case of a normal double-contact system pO| where e is proportional 
to the product ^T\r2- 

The average Josephson current-phase relation /(a) is shown to be strongly non-sinusoidal and roughly similar to 
the one known for a single nearly ballistic QPC, in the sense that it contains sharp "switching" between positive 
and negative values of the current as the phase varies via a = tt. The new feature of our system is that it is 
possible to vary the width of the swithching region ba by the electric gate potential Vg\ in the case of equal reflection 
probabilities r\ — this electric modulation is especially pronounced, ba oc | cos(7rCVJj/2e)|. The noise spectrum of 
the supercurrent is found to consist mainly of two peaks: the "zero-frequency" peak due to rare thermal exitations 
of the upper level of the system, and another one centered around the energy difference Wq. between the two levels. 
The widths of these peaks are determined by the inverse life-time r of the two states of our TLS, which is due to 
electron-phonon and electromagnetic couplings. Both these sources of level decay are expected to be very weak in 
the system considered, but the corresponding quantitative analysis is postponed for the future studies, so we present 
here only the results for the frequency-integrated (over those narrow intervals 1/t) noise power. 

The S-S-S device with almost ballistic contacts is a new type of a system which may be used as a realization of an 
artiflcial "spin 1/2" — an elementary unit for quantum computations. In comparison with usual Josephson systems 
with tunnel junctions which were proposed for the use in adiabatic quantum computatons [ll| ], the advantage of our 
system is that it may operate at considerably higher values of the Josephson critical currents; moreover, the current- 
phase characteristics of such a system is almost universal in the sense that it is determined mainly by the microscopic 
parameters of the SC materials and only weakly depends on the specifics of contact fabrication. 

We are grateful to K. A. Matveev, Yu. V. Nazarov and especially to G. B. Lesovik for many useful discussions. This 
research of M.V.F. was supported by the INTAS-RFBR grant # 95-0302, the coUaboration grant # 7SUP J048531 
from the Swiss National Science Foundation and the DGA grant # 94-1189. 
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